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ABSTRACT 

Observations of far-infrared rotational emission lines which arise in the 
shocked gas associated with Orion-KL are presented, including detections of the CO 
J = 34 -* 33, J = 31 -♦ 30, J = 26 -* 25, and OH 2 n 3 / 2 J p = 7/2” -* 5/2 + emission 
lines, as well as improved measurements of the CO J = 22 -» 21 and OH 2 n 3/2 J = 
5/2 -» 3/2 lines. These lines are observed to have velocity widths of Ai/ ~ 20-30 
kms -1 , somewhat less than either the 2 pm H 2 lines or the high-velocity "plateau” 
component of the millimeter-wave CO lines seen in this object. An H 2 column density 
of ~ 3 x 10 21 cm -2 , a total mass of ~ 1 ^ and characteristic temperature and den- 
sity T ~ 750 K and n(H 2 ) ~ 2 x 10 8 cm -3 can be derived from the CO intensities. 
The density is too low by at least an order of magnitude for the observed infrared H 2 
and far-infrared CO emission to be accounted for by a purely hydrodynamic shock, 
and support is lent to hydromagnetic shock models. 
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From the present measurements, the relative abundance of CO Is estimated to 
be C0/H 2 w 1.2 x 10 -4 , corresponding to 20% of the cosmic abundance of C existing 
in the form of CO. The average relative abundance of OH in the shocked gas is 0H/H 2 
£ 5 x 10 -7 , more than an order of magnitude greater than that expected in the 
quiescent molecular material surrounding the shocked region. The CO lines are opti- 
cally thin, and the OH lines have r £ 3. An upper limit to the intensity of the HD 
J = 1 -*0 line is used to derive an upper limit of 4x10 -5 for the D/H relative abun- 
dance in the Onon cloud core. 
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I. The Molecular Outflow in Orion-KL 

Orion-KL, the molecular region associated with the Kleinmann-Low nebula, is 
the prototype of the class of "active" molecular cloud cores. This recently- 
recognized class is characterized by hypersonic mass outflow, affecting large por- 
tions of the cloud core, which is apparently driven by mass loss from a newly-formed 
star. The young star itself Is generally manifested by a compact Infrared source, and 
the outflow displays three main features: motions of compact objects, molecular lines 
involving states of very high excitation, and molecular fines with very large velocity 
widths. First noticed were the motions of H 2 0 masers (e.g. Knowles eta/. 1969, 
Strelnitskii and Syunyaev 1973; Genzel and Downes 1977; Genzel etaJ. 1981) and 
of Herbig-Haro objects (e.g. Cudworth and Herbig 1979), which indicate large-scale 
expansion. Hot, shocked molecular matter associated with the outflows (7 £ 2000 
K) is seen in the near-infrared transitions of molecular hydrogen (e.g. Gautier et a/. 
1 976; Beckwith et al. 1 983). Molecular lines from these regions have widths as 
large as 100 km s -1 (e.g. Zuckerman, Kulper and Rodriguez-Kuiper 1976; Nadeau, 
Geballe and Neugebauer 1982); the spatial distribution of different velocity com- 
ponents of these lines often indicates a bipolar morphology (e.g. Snell, Loren and 
Plambeck 1 980; Bally and Lada 1 983). Large amounts of energy and momentum may 
be fed into the surrounding molecular clouds during the lifetime of the outflow, and 
these objects might, therefore, play a role in the support of molecular clouds against 
self-gravity (cf. Bally and Lada 1 983) and influence the rate of further star forma- 
tion in these clouds. 

To study the Interactions between the outflow and more quiescent material, 
we examine the hot shocked gas that lies at their Interface. The presently-available 
probes of this shocked gas include the rotation-vibration and pure rotational lines of 
H 2 in the near- and mid-infrared (cf. Beckwith et al. 1983, Beck 1984) and the far- 
infrared rotational lines of CO (Watson et al. 1980; Storey et al. 1981; Stacey et al. 
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1 982, 1 983) and OH (Storey, Watson and Townes 1 981 ). The far-infrared lines of 
CO are of particular importance because they are density-sensitive, and thereby 
carry information on the shock's compression ratio. 

In this paper we report observations of the far-infrared lines of CO and OH in 
the shocked matter associated with Orion-KL The present data include the first 
detections of CO J = 34 -+ 33 (X = 77.08 yum), J =31 -* 30 (X = 84.41 yum), CO J = 
26 -* 25 (X = 100.5 yum), and OH 2 n 3/2 J p = 7/2“ -* 5/2 + (X = 84.60 yum) lines, as 
well as improved measurements of the previously-observed lines CO J » 22 -» 21 (X = 
1 18.0 yum; Watson et a/. 1980) and OH 2 II 3/2 J = 5/2 -* 3/2 (X = 1 19.2, 1 19.4 /mi; 
Storey, Watson and Townes 1981). The CO measurements and the parameters 
derived from them are consistent with our previous observations (Watson et aJ. 
1 980, Storey et al. 1 981), but with the addition of the CO J = 34 -» 33 line provide 
stronger constraints on the density and compression ratio of the shock. The latter is 
shown to be too small for a simple, purely hydrodynamic shock to account for the 
observed activity. The data generally support the shock models of Chernoff, Hollen- 
bach and McKee (1982) and Dralne end Roberge (1982, 1984), in which magnetic 
effects are critically important to the shock's characteristics. 

II. Observations 

The instrument used In this experiment was the UC Berkeley tandem Fabry- 
Perot spectrometer (Storey, Watson and Townes 1980; Watson 1982). Observa- 
tions were made using velocity resolutions of 60 km s -1 (at 77 and 84 /im), 75 
kms -1 (at 100 and 112 yum), and 110 km s -1 (at 119 /tm); the instrumental profile 
is a Lorentzian. The system noise equivalent power (NEP) was 1-2 x 10 -14 W 
H z~ 1/2 , including all losses in the spectrometer, telescope and atmosphere. The 
observations took place in 1982 February and 1983 February on the 91.4 cm tele- 
scope of the NASA Kuiper Airborne Observatory. At the observing altitude of 1 2.5 km, 
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the water vapor column density along the line of sight was typically 20 microns. The 
system field of view was 44" in diameter (FWHM) and the solid angle of the beam 
was 5 x 10 -8 sr for all observations except those at 84 /mu, where the correspond- 
ing parameters were 30" and 2.7x1 0“® sr. The telescope's oscillating secondary 
mirror was used to provide a 4' east-west chop at 29 Hz. 

Each line was observed with the beam centered on the Becklin-Neugebauer 
object (a ig50 - 5 h 32 m 46. s 7, <5 ig50 = -5 24' 17"). Additional points in a 
northwest-southeast strip were observed in the 119 fjm OH lines. Line fluxes were 
derived from the measured line-to-continuum ratio and the photometric measurements 
of Werner etal. (1976). The accuracy of the fluxes thus obtained is about 30%; 
the accuracy with which the ratio of two line fluxes is determined is somewhat better 
than this (15-20%). At 84.5 /im, the different beamsize makes necessary a correc- 
tion for source coupling, which we estimated on the basis of our [O i] 63.2 fjm obser- 
vations in the direction of BN with both the 30" and 44" beams (see Werner et a/. 
1 984 and Ellis and Wemer 1 985). At this position, the [0 ij line is dominated by 
emission associated with the shocked region. The ratio of [O i] flux in the two beams 
is 1.5 (44" to 30"), and this factor has been applied to the 30" CO and OH measure- 
ments. 

The spectrometer's velocity scale was determined by observing fringes of 
0.6328 /m He-Ne laser light reflected from the scanning Fabry-Perot interferometer. 
Wavelength reference points were provided by nearby HDO, H 2 18 0 (Kyro 1981, 
Messer, Helminger and DeLucia 1 983), H Z S (Raud, Camy-Peyret and Johns 1 983) and 
NH 3 lines (Urban et aJ. 1981), which were measured using the spectrometer's internal 
absorption cell. The CO and OH wavelengths are obtained from measurements by 
Todd et al. (1 976) and Brown et a/. (1 982), respectively. 

Table 1 is a summary of the observations made with the beam centered on BN. 
The spectra of the individual CO and OH lines are shown in Rgures 1-5. The results 
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Table 1 

CO and OH Lines From Orion-KL 


Line 

X(jum) 

Flux 4 

(10“ 17 Wcm" 2 ) 

Intensity 41 * 

(1 0~ 2 erg s -1 cm -2 sr -1 ) 

CO J = 21 -*20 

124.194 

5.1 d 

1.0 d 


119.441 

0.7 ±0.3 

0.1 5±0.06 

0H 2 n 3/2 J = 5/2-*3/2 


119.234 

0.83 

0.17 

13 C0 J = 23-*22 

118.660 

<0.3 

<0.06 

CO J = 22 -* 21 

118.581 

5.3 

1,1 

CO J = 26 - 25 

100.461 

1.9 

0.38 

CO J = 27 -* 26 

96.7725 

3.2(+0.9,-1.9) 

0.64(+0. 18,-0.38)* 


84.5966 

0.56 c 

0.17° 

OH 2 n 3/ , J = 7/2-*5/2 


84.4202 

<0.2 C 

<0.05° 

CO J = 31 - 30 

84.4107 

0.6-0.87° 

0.18-0.26° 

CO J = 34 -* 33 

77.0586 

0.43 

0.086 


4 The uncertainty in the absolute fluxes and intensities is ±30% unless 
Indicated otherwise; the uncertainty in ratios of intensities masured 
with the same beam is somewhat smaller than this (15-20%). 

* Averaged over a 44" (FWHM) beam. 
c 30” beam. 

* Storey et al. (1981). 


of the additional OH observations are shown in Figure 6, along with similar "strip 
maps" in the lines of CO J = 21 -> 20 (Storey ef al. 1981) and H 2 v = 1 -» 0 S(1 ) 
(R.R.Treffers, personal communication; Beckwith et al. 1978). One of the spectra 
(Figure 2) also covers the J = 23->22 line of 13 CO; unfortunately, interference from 
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telluric 0 3 allows only a modest upper limit to be set for the Intensity of this line. 

The absolute frequency reference for the 84.5 / 4 m observations was obtained 
by assuming that the longer-wavelength line in Figure 4 is the OH 
2 n 3/2 J p = 7/ 2 + -»5/2 - line, and that its velocity centroid Is the same as the 119 
/ xm OH lines. In doing so, we ignore an apparent overall redshift of the spectrum sug- 
gested by a calibration with respect to the NH 3 J = 6-*5, s-»a lines produced in the 
spectrometer's internal gas cell. The apparent redshift would be 24 km s -1 (0.4 of a 
spectral resolution element), and could have been caused by a small misalignment of 
a mirror. At the spectral resolution of 60 km s -1 , the CO J =31 -»30 and OH 84.42 
/xm line are blended. However, the resulting line peaks near the CO wavelength, and 
the integrated intensity of this feature is not very different from that expected for 

the CO J = 31 -+30 line alone on the basis of the intensities of the other CO lines 
* 

observed. Thus * e consider the 84 /xm observations to have produced clear detec- 
tions only of the CO line and the 84.6 / 4 m OH line. We estimate that at least 70% of 
the intensity of the observed 84.4 /xm feature can be accounted for by the CO tran- 
sition, resulting in a limit on the OH 84.4 /x m/84.6 /xm ratio of <0.5. (The 84.4 /xm 
feature does seem to show some excess emission on the redshifted side; this may 
represent the contribution of the OH transition). 

At 1 1 9.4 /xm, the 2 n 3/2 J p - 5/ 2 + -*3/ 2~ line is not resolved from the tellu- 
ric 0 3 line at 1 19.46 /xm, and the OH line emission must be corrected for continuum 
emission absorbed by this line. Another telluric 0 3 line at 1 19.05 /xm (see Figure 1), 
which has 1 .9 times the equivalent width of the former 0 3 line, provides us with an 
estimate of this absorption. With the correction, similar intensities for the 1 1 9.4 and 
1 1 9.2 /xm OH lines are found. 

In addition to the CO and OH data, upper limits to the intensities of the HD 
J = 1 -*0 line and several lines in the J = 5-*4 and J = 6-»5 multiplets of NH 3 were 
obtained, and are listed in Table 2. The NH 3 data are consistent with the 
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J = 4-*3, a -*s, K = 3 detection reported previously (Townes eta/. 1983), and will 
not be discussed further here. The HD line was observed at 75 km s -1 resolution. 
Interference from telluric HDO, H 3 7 0 and 0 3 lines produced a baseline uncertainty of 
2% of the continuum at 1 1 2 pm, leading to the upper limit shown In Table 2. 


Table 2 

HD and NH 3 Upper Limits 



Une 

A 

(pm) 

Flux 

(10~ 17 Wcm -2 ) 


HD J = 1 -*0 

1 1 2.075 

<0.1 a 

NH 3 j 

= 5-»4, a -*s, K = 1 

100.277 

<0.5 a 

nh 3 j 

= 5-»4, a -*s, K = 2 

100.213 

<0.5 a 

NH 3 j 

= 5-*4, a-*s, K = 3 

100.105 

<0.5 a 

nh 3 j 

= 6-*5, a -*s, K = 4 

84.6198 

<0.2^ 

nh 3 j 

= 6-* 5, a-*s,K= 5 

84.5422 

<0.2* 


a 44" beam. 

b 


30" beam. 
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Wavelength (pm) 

1194 1192 1190 



+ 100 0-100 +100 0-100 


v L sr (km s' 1 ) 


Figure 1: OH 2 U 3/2 J = 5/2 -* 3/2 emission in Orion-KL. This spectrum is the result 
of 7 minutes of integration (6.6 s per point). As discussed in the text, the actual in- 
tensities of the two OH lines are very similar after accounting for the effect of the 
telluric 0 3 tine at 1 19.46 /jjm. The wing of a telluric 0 2 line centered at 1 19.81 /urn 
has been divided out of the spectrum. 

Wavelength (pm) 

1187 1186 1185 



v Lsn (kms* 1 ) 

Figure 2: The CO J = 22 -♦ 21 line In Orion-KL. The total integration time for this 
spectrum was 3 minutes, or 3 s per point. An upper limit to the intensity of the 
13 CO J = 23-* 22 line was also derived from this observation. 
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Wavelength (pm) 

Figure 3: Detection of the CO J = 26 -♦ 25 line in Orion-KL. Note that the line Is 
clearly broader than the unresolved telluric H 2 18 0 line. Upper limits to the intensities 
of three NH 3 lines were also obtained from this spectrum. The total Integration time 
was 30 minutes (14 s per point). 


Wavelength (pm) 


84 6 84 5 84 4 



Figure 4: Detection of the OH 2 n 3/ , 2 = 7/2“-*5/2 + and CO J = 31 -*30 lines In 

Orion-KL. The velocity scales are appropriate for the CO line and the 84.6 pm OH line, 
the latter having been taken as the absolute wavelength reference for the observa- 
tion. At this velocity resolution, the CO line is blended with the OH 2 n 3/2 
J p = 7/ 2 + -*5/2 - line, but, as discussed in the text, £70% of the observed line is 
probably due to the CO transition. Upper limits to the intensities of two NH 3 lines 
were also obtained from this spectrum. The total Integration time was 22 minutes 
(21 s per point). 
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Wavelength (pm) 



Vlsr (km s' 1 ) 

Figure 5: Detection of the CO J = 34 -* 33 line in Orion-KL. The total 
Integration time for this spectrum was 42 minutes (40 s per point). 
Baseline curvature, caused by strong absorption lines centered on ei- 
ther side of the spectral interval shown above, was removed by divid- 
ing by the observed spectrum of Mars. 

III. Physical Parameters and Abundance of CO in the Shocked Gas 

a) Velocity Distribution 

The CO line observations were done at high enough spectral resolution to par- 
tially resolve the lines (see Figure 3 in particular). The observed linewidths are all 
1 5-25 km s -1 broader than the instrumental full width at half-maximum. This prob- 
ably indicates that the far-infrared CO lines are narrower than the plateau component 
of the millimeter-wavelength lines of CO (Ai/ a 50 km s'* 1 ; e.g. Knapp et al. 1981) or 
the H 2 lines at 2 fjm (Ai/ a 40 km s -1 ; Nadeau, Geballe and Neugebauer 1982, Sco- 
ville et al. 1 982). To show this, we have performed numerical convolutions of these 
latter lineshapes and our instrumental profile. If the CO plateau profile, observed at 
very high velocity resolution with beams similar to ours, is taken, the resulting far- 
infrared CO profile is £ 50 km s -1 wider than the instrumental FWHM. In the case of 
the 2 yum H 2 lines, we take the ‘'intrinsic" H 2 v = 1-*0S(1) profile presented by 
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0 V 2' 3' 4" 



Figure 6: Strip maps of H 2 , CO, and OH emission in Orion-KL. The strip 
lies along a line leading 30 degrees W of N through H 2 Pk 2 (Beckwith 
eta/. 1978). The H 2 strips are made with 45" resolution (bold curve) 
and 1 3" resolution (thin curve), and come from Treffers (personal com- 
munication) and Beckwith etal. (1978), respectively. Also shown on 
the H 2 strip map are the positions of Pk 1 and Pk 2. 


Nadeau, Geballe and Neugebauer (1982). This profile results from a weighted aver- 
age of data from many small-beam measurements covenng the bulk of the shocked 
region, from which the effect of the finite spectrometer resolution (20 km s _1 ) was 
removed by numerical deconvolution; it still has the effects of extinction within the 
shocked region, which makes the observed line narrower than the true velocity 
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distribution. When convolved with our instrumental profile, the "intrinsic" H 2 profile 
gives rise to far-infrared CO linewidths which exceed the resolution by 
40-50 km s -1 , substantially greater than those observed. Since the far-infrared CO 
lines are not affected by extinction, they would appear still wider than this if they 
resulted from the same velocity distribution as the lines. Thus the sources 

of the CO plateau and the 2 /xm H 2 lines probably do not contribute much of the 
observed far-infrared CO emission. 

If the far-infrared CO lines have the same shape as the millimeter-wave CO 
lines, but a different width, the intrinsic width lies in the range 20-30 km s -1 
(FWHM). 

The OH lines are also observed to be about 1 5-20 km s -1 broader than the 
resolution, and thus may have an intrinsic velocity dispersion similar to that of the 
far-infrared CO lines. 

b) Density, Temperature and CO Abundance 

The upper states of the nine far-infrared CO lines observed to date in Orion- 
KL range in excitation temperature Ej/k from 750 K (J = 1 6) to 3300 K (J=34). Thus 
the manner in which the shocked gas is sampled by these lines is similar to that for 
the H 2 v = 0 S(2) line at 1 2.28 pm, for which E d /k = 1 700 K. In contrast to the H 2 
rotational ladder, which is populated according to thermal equilibrium, the far-infrared 
CO observations involve states which are rather far from thermal equilibrium as well 
as some which are fairly close, so the CO line ratios are sensitive to density as well 
as temperature in the shocked gas (Storey et al. 1981). Combination of thi3 informa- 
tion and that obtained from observations of the H 2 v = 0 S(2) line yields the C0/H 2 
relative abundance and mass of the shocked gas. 

The method by which the H 2 density, temperature and CO/H 2 relative abun- 
dance are estimated is discussed by Storey et a/. (1981); we outline it briefly here. 
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Sin ce the far-infrared CO lines in Orion-KL are unextinguished and optically thin (cf. 
Watson et al. 1 980), their intensities are expressed by: 




hv 

47T 


AjNj, 


( 1 ) 


where A j is the A -coefficient for the J -* J-1 transition and Nj is the column density 
of CO in the J th state. Given the molecular hydrogen density n(H 2 ) and temperature 
T, the Nj's are determined by: 


N, 


Aj + n(H 2 ) ^7 jj- - N J+ i 4 J+1 + n(H 2 ) Nj. y j.j 


( 2 ) 


Y j N J =N { 


CO ' 


where jjj> is the rate coefficient for the collisional transition J -* J' and N co is the 
total CO column density. It is assumed here that decay of higher vibrational levels 
pumped by collisions or by radiation in the 4.6 pm CO fundamental vibration bands has 
a negligible effect on the populations of the rotational levels. This assumption may be 
made plausible by considering the vibrational collisional excitation rates (Scoville, 
Krotkov and Wang 1 980) and number of available 4.6 /xm photons (cf. Downes et al. 
1981). Values of n(H 2 ) and T are obtained by finding the solutions to Equations 2 
which match the shape of the Nj distribution inferred from the observed fluxes. Once 
n(H 2 ) and T are known, Nj/N co is determined for all J , and N co may be obtained from 
any of the observed line intensities by using Equation 1 . For the molecular constants 
in Equations 2, we take the collisional rate coefficients jjj. calculated by S. Green 
(McKee et al. 1982), and compute the A -coefficients using the rigid-rotor approxi- 
mation and a permanent dipole moment of 0.1 1 2 D . 1 

® Green and Chapman (1983) and Chackerlan and Tipping (1983) have computed higher-order 
corrections to the collisional rate coefficients and the A -coefficients (respectively). These 
corrections are not Included In the results presented here because they are relatively small and 
nearly cancelling Upon repeating our calculations using these corrections, we found the only ef- 
fect to be a 30% reduction In the derived H 2 density. 
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Much is already known about the Orion-KL shocked gas from H 2 observations. 
The latest observations can be summarized as follows. Rotation-vibration transitions 
at 2 /*m and highly-excited pure rotational transitions arise from matter at about 
2000 K; this gas lies behind an extinction of 2 magnitudes at 2 jum (Beckwith et al. 
1983) and its column density, averaged over a 1' diameter area centered on IRS2, is 
2.4 x 1 0 19 cm -2 . The H 2 v = 0 S(2) line at 1 2.28 /zm has also been observed in this 
region (Beck, Lacy and Geballe 1979; Beck eta/. 1982; Beck 1984). A value for 
the extinction at 12 /Mn of 0.75 mag can be derived from the 2 /an extinction using 
standard interstellar extinction curves (Becklin eta/. 1978). From this, and the 
observed intensities of the 12.28 /an line, we infer a column density in the v ~ 0, J = 
4 state of 1.5 x 10 20 cm -2 , averaged over the same 1' diameter area. This obvi- 
ously cannot be due to the 2000 K gas alone and indicates the presence of an addi- 
tional amount of somewhat cooler gas. To account for these two batches of gas in 
the interpretation of the far-infrared CO data, Storey et al. (1981) modelled the 
region with two components, one at 2000 K (designated "hot”) and the other one 
cooler (designated "warm"). The H 2 column densities derived from the 2 /un and 12 
/xm observations were used along with the solution of Equation 2 and the assumption 
of pressure equilibrium (so that r?(H 2 ) a 1/7) to fit the CO relative intensities and 
thus derive the temperature of the cooler component, the CO/H 2 relative abundance, 
and the molecular hydrogen density. The present CO observations are more exten- 
sive than those available at the time of the earlier calculation, and the accepted 
value for the average extinction at 2 /m has since changed from 4 mag to 2 mag. 
Therefore we have repeated the calculation. The range of parameters which produce 
an acceptable fit is shown in Table 3; the best fit is given by 
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7 = 7 50 K 

n(H 2 ) = 3x10® cm -3 
W(H 2 ) = 3x1 0 21 cm" 2 
CO/H 2 = 1.2X10" 4 

for the warm component. The CO intensities resulting from the "best fit" parameters 
is displayed in Figure 7. Note that the "hot" component, which dominates the 2 / 2 m H 2 
emission, contributes very little to the observed far-infrared CO lines, so that ours is 
essentially a single- component model of the "warm" post-shock gas. 



Figure 7: Observed CO line intensities, compared to the simple model 
discussed in the text. The intensity distribution for T = 750 K, 
A/(H 2 ) = 2.8x1 0 21 cm -2 , and thermal equilibrium (LTE) is also shown. 
The intensity plotted for the J - 31 -»30 line (square symbol) is the full 
intensity of the 84.4 yum feature, part of which may be due to the 
nearby OH transition. The 1 7 -» 1 6 and 1 6 -» 1 5 points are from Sta- 
cey et at. (1982, 1983). 
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Table 3 

Two-component Model: 
Orion-KL Post-shock Gas 

Parameter 

"Hot" gas "Warm" gas 

Temperature 

2000 K 1000-500 K 

n(H 2 ) (cm -3 ) 

0.5-3x10® 1-16x10® 

A/(H 2 ) (cm -2 ) 

2.4 xIO 19 2 - 6 x 1 0 21 

M( H 2 ) (1^ 

0.01 1 - 3 

co/h 2 

1.5 - 0.8 x 10 -4 1.5 - 0.8 x 10 -4 


The more important results of the far-infrared CO observations and the simple 
model described above are the limits placed on the CO/H 2 relative abundance, the H 2 
column density and total mass of the post-shock material, and the approximate H 2 
density. A single-temperature model for the warm component is of course an oversim- 
plification, since there must be a temperature gradient created by the cooling of the 
post-shock gas. However, we also note that the derived values of most of the 
parameters do not change radically over the temperature range where the fit is good, 
which indicates that similar values would be derived If the source Is considered to 
have temperature variations in this range. The post-shock molecular hydrogen den- 
sity cannot be much higher than a few times 1 0 6 cm -3 . The CO/H 2 ratio is a more 
accurate estimate, regardless of the precise values of density and temperature, 
because it comes from the flux ratio of two lines (CO J = 22 -» 21 and H 2 v = 0 S(2)) 
which are optically thin and arise from states which are nearly thermalized and have 


similar energies. 
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IV. Thermal OH Emission from the Shocked Gas 

The OH lines reported here are the lowest-lying pure rotational transitions in 
the 2 n 3/2 electronic state (see Figure 8). The 119 p,m OH observations (Figures 1 
and 6) confirm the tentative detection by Storey, Watson and Townes (1981). From 
the strip map of the 1 1 9.2 pm line (Figure 6) it appears that the source of the emis- 
sion Is about 1 1 in size (assuming a Gaussian source) and covers roughly the same 
region as the far-infrared CO and near-infrared H 2 emission. With the 119.4 pm 
Intensity corrected as discussed above (§ II), no variation of the 119.4/119.2 pm 
intensity ratio with position is seen. At the position of BN, the intensities of the four 
OH lines observed are in the proportions 0.9 : 1.0 : 1.0 : <0.5 (1 19.4 : 1 19.2 : 84.6 : 
84.4 p m). 


Energy (cm"’) 
300 


200 


J p 
L 7/2 : 


ioo h 


T 


84 59665 pm 


5/2 Z 


84 42015 pm 


119 4410 pm 


L 3/? : 


119 2336 pm 

L_ 


z n , /2 

p j 
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Figure 8: The lowest rotational energy levels of OH, labelled by rota- 
tional quantum number J and parity p, showing the observed transi- 
tions. The A-type splitting is exaggerated in this drawing. 


As we will see below, coilisional excitation and radiative trapping in the 
shocked region can reasonably account for the OH emission; purely radiative mechan- 
isms cannot. Both far-infrared and near-infrared pumping can be ruled out by consid- 
ering the luminosity of the continuum sources and the branching ratios of the 
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resulting cascades; the continuum at 53 and 37 /m is too weak by at least an order 
of magnitude to produce the observed lines through a rotational cascade, while the 3 
fjm continuum sources are at least four orders of magnitude too weak to do so with a 
vibrational cascade. (Of course, direct rotational pumping by 1 1 9 and 84 fjm photons 
cannot work because the OH lines would then appear in absorption rather than emis- 
sion.) Resonant scattering of far-infrared continuum photons by optically thick OH 
lines also cannot produce the observed spectra; given the color temperature of dust 
in the Orion-KL region (80-100 K), resonantly-scattered continuum radiation would 
produce lines with relative intensities 1 : 1 : £2 : ^2, contrary to observations. 

There are considerable uncertainties in the collisional rate coefficients for OH 
rotational excitation, especially in the relative rates for different A-doubiing levels 

m 

with the same angular momentum quantum number J (Andresen ef a/. 1 984a, b), so 
any numerical results derived from the OH intensities must be regarded as very 
approximate. Nevertheless, it seems clear that the OH relative abundance in Orion-KL 
is unusually large. The relative abundance of OH may be estimated from the 1 1 9 /^m 
line fluxes as follows. Since these lines have rather large A -coefficients (0.138 
s -1 ) and involve the ground rotational state, one might expect them to have large 
optical depth but small fractional upper-state population. Though the former condi- 
tion would imply that the photons are absorbed and re-emitted many times on their 
way out of the cloud, the latter means that collisional de-excitations are not fre- 
quent enough to destroy a significant number of them. Thus the line flux should be 
proportional to the rate of collisional excitations (or, equivalently, the molecular emis- 
sion measure yVj 0H n(H 2 )ds), just as it would be if the line were optically thin and 
subthermally excited, though there would be an additional geometry-dependent fac- 
tor due to the repeated scattering of the line photons. Assuming plane-parallel 
geometry and the validity of the escape probability approximation, the line intensity 
lj is given by (see Appendix A): 
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>J - W 0H (3) 

From the observed strength of the pair of 119 /a OH lines, the density of 
rr(H 2 ) = 3x10® cm -3 derived In the previous section, and the theoretical collisional 
excitation rate coefficients at T - 600 K (Dewangan and Flower 1 982), we obtain 
the following values of OH column density and line optical depth averaged over the 
beam: 


< W 0H > = 1 .5x1 0 15 cm -2 , 

< r > = 3 (each line) . 

A line width of 30 km s -1 was taken in the calculation of the average optical depth. 

«■ 

Both of these average values should be regarded as lower limits to the actual quanti- 
ties because of the assumption of plane-parallel geometry and since the source is 
not likely to uniformly fill the beam. As may be easily shown, these results imply small 
excited-state populations as long as the beam filling factor is larger than about 0.05. 

The OH emission appears somewhat more localized around Pk 1 than does the 
CO; the two may not be identically distributed. However, if the distributions are 
taken to be approximately the same, a relative OH abundance of 0H/H 2 = 5 x 10 -7 is 
implied. Because of the assumptions of plane-parallel geometry and distribution 
throughout the shocked layer, this is a lower limit to the local OH/H 2 ratio, yet it is 
more than an order of magnitude greater than the 0H/H 2 ratio under any conditions 
which are likely to be found in dense, quiescent molecular clouds (Prasad and Hun- 
tress 1 980). Such an overabundance is the expected outcome of the elevated tem- 
perature in a post-shock molecular environment, in which endothermic or high- 
threshold reactions between neutral molecules can occur; in particular, the reaction 
O + H 2 -♦ OH + H (Hollenbach and McKee 1 979, 1 980). 
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V. Limits on the D/H Ratio in OMC-1 

The iine-to-continuum ratio of the J = 1 -*0 transition of HD can be used to 
obtain the HD/H 2 abundance ratio (= 2D/H) in dense molecular clouds such as Orion. 
Because HD Is expected to bear most of the deuterium in molecular clouds, this 
measurement of D/H is free of the effects of chemical isotopic fractionation which 
introduce uncertainties in D/H determinations by millimeter-wavelength lines of trace 
species. Since the J =1 state is only 1 28 K above ground, emission from the direc- 
tion of Orion-KL would be dominated by the warm, quiescent core of OMC-1, for which 
T gas ~ Tdutt ~ 70 K, A/(H 2 ) ~ 10 24 cm -2 and r D ~ 0.3 for the 100 fjm continuum emis- 
sion from the central 1' (cf. Werner 1982). In Appendix B, the following expression 
is derived for the HD/H 2 ratio in such clouds: 


HD 

Ai/ c 

1 11 7, 22 K 

1 

ls T ° - "i] 

l( 1-0) 

h 2 

V 

13K 7 T 

•J _ q— 128.4 K / T 

! T D 



where hv c /v is the effective fractional continuum bandwidth of the spectrometer 
(.nhv FWHM /2v for a Fabry-Perot interferometer) and /(I -*0) is the integrated inten- 
sity of the line. With the above values of temperature and continuum optical depth, 
the upper limit of 0.02 for the line-to-continuum ratio implies HD/H 2 < 8x1 0 -6 , or: 

-5- <4x10" 5 . 
n 

The interstellar abundance of deuterium, which is thought to be near 10~ 5 , is not 
very strongly constrained by this measurement. 
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VI. Comparison with Theoretical Models 
Although the density and temperature estimates obtained from the far- 
infrared CO observations are not very precise, they are good enough to rule out the 
possibility that a purely hydrodynamic (non-dissociative non-magnetic, or NDNM) 
shock is responsible for the H 2 and CO emission in Orion-KL, as may be seen quite 
simply. Such a shock would have an initial compression of a factor of 6, with an 
accompanying temperature jump from T 0 to T m ax . The post-shock gas would have 
approximately constant pressure, so that the post-shock density would be n = 6 n 0 
Tmax/f. Near-infrared observations of H 2 lines (e.g. Beckwith eta/. 1983, Beck and 
Beckwith 1983) and conservation of energy provide constraints on the maximum 
temperature and preshock density: 


7" max 2700 K , 

n 0 Z 47r/(H 2 )/m(H 2 ) v 3 £ 10 a cm~ 3 , 

since the shock speed v s must be <; 15 km s'" 1 to keep the molecular hydrogen from 
being collisionally dissociated (cf. Kwan 1977). By the time the post-shock gas 
would have cooled to 750 K, the density would have been £ 2x1 0 7 cm -3 , and thus 
all of the states from which the far-infrared CO lines onginate would be thermalized 
at high temperature, in conflict with the observations. 

Because of this difficulty with NDNM shock models, and the additional problem 
posed by the very large velocity widths seen in the near-infrared H 2 lines (Nadeau, 
Geballe, and Neugebauer 1 982, Scoville et at. 1 982), theoretical attention has 
turned to hydromagnetic shocks (Draine 1980, Draine, Roberge and Dalgarno 1983, 
McKee, Chernoff and Hollenbach 1984, Draine and Roberge 1984). These shock 
structures have been applied specifically to Orion-KL by Chernoff, Hollenbach and 
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Table 4 

Comparison of Models and Observations 
Wavelength Observed intensity Calculated intensity 6 


Line 


(10 -2 erg s -1 cm : 


(/xm) 

(10 -2 erg s -1 cm -2 sr -1 ) 



CHM DR 


CO 

124.2 

1.0 a 

0.90 

0.98 

OH 

119.4 

0.15 

0.18 

0.17 

OH 

119.2 

0.17 

0.18 

0.17 

CO 

118.6 

1.1 

... 

0.92 

CO 

100.6 

0.38 

... 

0.56 

CO 

96.77 

0.64 a 

0.40 

0.55 

CO 

87.19 

0.24 a 

0.22 

0.22 

OH 

84.60 

0.17 

0.22 

... 

OH 

84.42 

<0.08 

0.22 

... 

CO 

84.41 

0.18-0.26 

... 

0.15 

CO 

77.06 

0.086 

0.10 

0.068 


* Storey eta/. (1981). 

6 CHM s Chernoff, Hollenbach and McKee (1982); DR a Drains and 
Roberge (1982). 

McKee (1982) and Draine and Roberge (1982). Calculated line intensities from 
these models are compared with the present far-infrared observations in Table 4, and 
the model parameters are listed in Table 5. In both cases, close agreement with 
experiment is obtained with "C-type" shocks (no discontinuous "jumps" in the fluid 
parameters). 
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Table 5 

Model Parameters 


Parameter 

CHM 

DR 

n 0 (H 2 ) (cm -3 ) 

1 x 10 5 

3.5 x 1 0 5 

^Zfun (mag) 

2.5 

4 a 

A 12 pm (mag) 

0.94 

0.8 a 

i/ s (km s -1 ) 

36 

38 

co/h 2 

6x1 0 -4 

1 .4 x 1 0 -4a 

B 0 (mG) 

0.45 

1.5 


a Assumed in all model calculations. Other values result from tree- 
parameter fits. 


VII. Summary 

The more important results of the far-infrared molecular line observations of 
Orion-KL are the following. 

1. Downstream from the 2000 K shocked gas seen in the 2 /zm H 2 lines, there is a 
much larger, cooler component of post-shock gas with column density A/(H 2 ) * 3 x 
10 21 cm -2 , total mass M * 1 Mq and characteristic temperature and density T * 750 
K and n(H 2 ) a 3 x 10 6 cm -3 . The density is too low by at least an order of magni- 
tude for the shock to be non-dissociative and non-magnetic. 
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2. The far-infrared CO and OH lines arising in this shocked gas have velocity disper- 
sions of Av « 20-30 km s _1 (FWHM), somewhat less than either the 2 /xm H 2 lines or 
the high-velocity plateau component of the millimeter-wave CO lines. The far- 
infrared CO lines are optically thin; the 1 1 9 pn OH lines have r £ 3. 

3. The relative abundance of CO is C0/H 2 =» 1.2 x 10 -4 , corresponding to about 20% 
of the cosmic abundance of C tied up in CO. Since no significant creation or destruc- 
tion of CO is expected to take place in the shock, this probably represents the nor- 
mal CO relative abundance in the Orion cloud, and can be compared with the value of 
£ 12% determined by Dickman (1978) from a 13 C0 dark-cloud survey. 

4. The average relative abundance of OH in the warm postshock gas is 0H/H 2 £ 
5x1 0 -7 , more than an order of magnitude greater than that expected in the quies- 
cent pre-shock gas. 

5. The relative abundance of deuterium in the core of the Orion molecular cloud is 
D/H < 4x10 -5 . 

We wish to thank the staff of the Kuiper Airborne Observatory for their excel- 
lent support of this experiment. We are also grateful to S.C. Beck, D.F. Chernoff, B.T. 
Draine, D.R. Flower, S. Green, DJ. Hollenbach, C.F. McKee, E. Serabyn and R.R. 
Treffers for discussing their results with us, and to M.K. Crawford for assistance with 
the observations. This project was supported in part by NASA grant NGR 05-0p3- 
511. 

Appendix A 

Derivation of Equation 3 

Let us assume the validity of the escape probability approximation and con- 
sider the transition k -*J in a molecule for which the absorption coefficient of the line 
is much greater than that of the dust. (This includes the lower rotational transitions 
of the most common molecules except those of molecules with zero or extremely 
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small dipole moments, like H 2 or HD). The intensity l kj emitted in this line, in excess 
of the continuum, is given by 

l kJ = fdv<p(v)\s kj (1 -p kJ )-B v <J D )(1 — e~^ D )] , (A1) 

where is the line profile and and (3 kj are the source function and escape 
probability for the line. We can neglect the term due to dust emission and absorption, 
since in the present application (shocked gas), the dust optical depth t d is small and 
the dust temperature T 0 is much lower than the gas temperature. In the limit of small 
population in the upper state k, the source function is simply 




2h v 3 



n j 9k 


(A2) 


where the g 's are the degeneracies of the given states. In the same limit, with state 
J as the ground state, the equations of statistical equilibrium which determine the 
population of state k reduce to 


n j n X 7 J‘ = n K- n(H 2> 7kj + n k Pkj A kj • (A3) 

/ ak 


Under the assumption that the effective radiation rate fi kj A kj is much larger than the 
collisional de-excitation rate n(H 2 ) y kJ , we obtain 


l kJ - f dvpiv) 


2ft i/ 3 9 ± 
c 2 9k 


] «( H 2 ) 7ji 

:k 

Pkj A kj 


(A4) 


In a plane-parallel cloud with a large velocity gradient, the escape probability is as 
given by Scoville and Solomon (1974): 


Pkj 


1 -e- 3T * 

3 Tjk 


sa 


1 

3 T jk 


(A5) 
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where the optical depth r jk (assumed large) is given by 


f dv<p{v) T Jk = ^ ( Nj B Jk - N k B kJ ) 


" T N » B * 


(A 6 ) 


- — ~ ", 


Bin? 9j 


(Note that under the approximations we have made so far, A/ totaj w Nj). Rnally, tak- 
ing the line profile 55 ( 1 /) to be narrow, we obtain 


't, = <•(«,)«, ■ 


(A7) 


This expression is identical to Equation 3, which is used above to compute the OH 
column density from the ground-state rotational lines at 1 1 9 fim. 

Appendix B 

Derivation of Equation 4 

The HD molecule has a very small permanent dipole moment {ji = 5.9x1 0 -4 D; 
Trefler and Gush 1968). Thus its rotational energy levels reach thermal equilibrium at 
very low densities, and its rotational lines are optically thin. In fact, the transfer of 
HD line radiation in molecular clouds is controlled by dust absorption, since even for 
very narrow lines the HD line absorption coefficient, given by 


K J -J -i 


8 rrW-i,j v 9 j L 
3hc Bv gj_^ nj -' l 


-hw/kT 

> 


(B1 ) 


is much smaller than that of interstellar dust, for which 


K 0 W 


10 


24 

n(H 2 ) ' 

125 fm 


1 cm -3 , 

A 


(B2) 
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(Savage and Mathis 1979, Whitcomb et al. 1981). Giant molecular clouds often have 
substantial far-infrared continuum optical depth, and this must be accounted for in 
the interpretation of HD intensities. Here, we will once again assume the validity of 
the escape probability formalism, and compute the HD J = 1 -*0 line intensity. Con- 
sider a dense, quiescent cloud in which the gas and dust are in thermal equilibrium, so 
that the source function is equal to the same Planck function for the HD line and 
the dust continuum. The escape probability for continuum photons is approximately 

/S D = e where r D = - fields is the continuum optical depth, and the escape pro- 
bability for line photons can be written as where (3 L is the probability for a pho- 
ton to escape line absorption. Since the self-absorption of HD is very small, we take 
ft = 1 - r L , where t l - - f^^ds is the portion of the total optical depth due to 
the line. Then the excess of specific line intensity above the continuum can be writ- 
ten as 

A/, = e u (n(i -ft/? D ) -s„(n(i -/? 0 ) 

= T L B v (T)e~ r ° . (B3) 

The present measurements would not spectrally resolve lines originating in quiescent 
clouds, so we compute the integrated line intensity /(I -*0): 

/(I -«0) = Jdvf{v)ZJ v - e~^°fd v fds B u k^ 0 

= e ~ T °fd v fds j i _ 0 (in LTE) 

= ^i-o'iW H oe‘ r ° . (B4) 

In thermal equilibrium, the fractional population of the J = 1 rotational state is ade- 
quately approximated by 
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3e- hu/kr 

f '~ 2W- 1 h v 

hv 3 30kT 

The intensity of dust continuum adjacent to the line is given by 

1 -e^Ai/c , 


(B5) 


(B6) 


where Au c is the continuum bandwidth of the spectrometer. By combining Equations 
B2, B4, B5 and B6, we obtain 


HD _ M T , 22 K ) [ 1 fe T ° — l](/(1 -»0) 

- [ ~ J [?3k +1 ' 7+ — J (i - e -iSEfTGTj J [- 7 ^ 


(B7) 


This expression is used above (section V) to derive an upper limit for the D/H relative 
abundance in core of OMC-1. 
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16 Abstract 

Observations of far-infrared rotational emission lines which arise m the shocked gas asso- 
ciated with Onon-KL are presented, including detections of the CO J = 34 -*• 33, J = 31 -*• 30, 

J = 26 -*• 25, and 0H 2 n 3 / 2 J 9 = 7/2” ■+ 5/2 emission lines, as well as improved measurements of the 

CO J = 22 -*■ 21 and OH 2 n 3 / 2 J = 5/2 -*• 3/2 lines. These lines are observed to have velocity 

widths of Av - 20-30 km s -1 , somewhat less than either the 2 ym H 2 lines or the high-velocity 
"plateau" component of the millimeter-wave CO lines seen in this object. An H 2 column density of 
~3 x 10 21 cm -2 , a total mass of -1 M e and characteristic temperature and density T * 750 K and 
'2 x 10 6 cm -3 can be derived from the CO intensities. The density is too low by at least an order 
of magnitude for the observed infrared H 2 and far-infrared CO emission to be accounted for by a 
purely hydrodynamic shock, and support is lent to hydromagnetic shock models. 

From the present measurements, the relative abundance of CO is estimated to be 
C0/H 2 1 1.2 x io 4 , corresponding to 20% of the cosmic abundance of C existing in the form of CO. 
The average relative abundance of OH in the shocked gas is 0H/H 2 > 5 x 10~ 7 , more than an order of 
magnitude greater than that expected in the quiescent molecular material surrounding the shocked 
region. The CO lines are optically thin, and the OH lines have x > 3. An upper limit to the 

intensity of the HD J = 1-K) line is used to derive an upper limit of 4 x 10 -5 for the D/H rela- 

tive abundance in the Orion cloud core. 
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